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ABSTRACT 



We analyze the impact of Fermi gamma-ray observations (primarily non- 
detections) of selected nearby galaxies, including dwarf spheroidals, and of clus- 
ters of galaxies on decaying dark matter models. We show that the fact that 
galaxy clusters do not shine in gamma rays puts the most stringent limits avail- 
able to-date on the lifetime of dark matter particles for a wide range of par- 
ticle masses and decay final states. In particular, our results put strong con- 
straints on the possibility of ascribing to decaying dark matter both the increas- 
ing positron fraction reported by PAMELA and the high-energy feature in the 
electron-positron spectrum measured by Fermi. Observations of nearby dwarf 
galaxies and of the Andromeda Galaxy (M31) do not provide as strong limits as 
those from galaxy clusters, while still improving on previous constraints in some 
cases. 



1. INTRODUCTION 



A large variety of independent astronomical observations have revealed that most of 
the mass in the universe is some form of non-baryonic, non-luminous, cold dark matter of 
as yet unknown composition. Understanding the fundamental nature of dark matter is one 
of the biggest o utstanding pro blems in cosmology and particle physics (for a comprehensive 



review, see e.g. iBertond l2010h . Compelling particle candidates for the dark matter, collec- 
tively termed weakly interacting massive particles (WIMPs), exist in several well-motivated 
theoretical extensions of the Standard Model of particle physics like s upersymmetry (light- 
est supersymmetric particle, for a review see e.g. lJungman et al.lll996) and Universal Extr a 
Dimensions (lightest Kaluza-Klein particle, for a review see e.g. iHooper fc Profumo! 120071 ). 
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A signal from dark matter could be detected in astronomical observations from WIMP pair 
annihilation or decay. For example, WIMP annihilation or decay generically leads to the pro- 
duction of gamma-rays (as w ell as of other high e nergy particles), potentially at a detectable 
level with current telescopes (IBertone et al.ll2005l ). The unparalleled sensitivity of the Fermi 
Large Area Telescope (LA T) at GeV energi es makes it an excellent instrument to look for 
this gamma-ray signal (e.g. lBaltz et al.ll2008l ). In addition, the recent detectio n by PAMELA 



of a p ositron excess relative to cosmic-ray secondary production in the Galaxy ( lAdriani et al. 
2009 1 ) and a gentle but clear feature in the Fermi-LAT positron-electron (e + — e~) spectrum 



(lAbdo et al.l 120091 ) has led to intense discussion of the possibility that one or both of these 
signals are due to dark matter. If true, we expect a corresponding signal in gamma-rays 
stemming from both the annihilation or decay event itself, as well as from the secondary 
inverse Compton (IC) up-scattering of background radiation by the high energy e + and e~ 
produced. 

As far as dark matter annihilation is concerned, first year Fermi observations of dwarf 
spheroidal galaxies, clusters of galaxies, the extragalactic gamma-ray background, and searches 
for gamma-ray lines have placed some of the strongest constraints yet on the possible WIMP 
pair annihilation cross-section and have alrea dy excluded some annihilation models that 



might explain the PAMELA positron excess (lAbdo et al.l l2010bl ; lAckermann et al.l 12010b 



Abdo et al.ll2010dl ial). In addition, several authors have recently pointed out that gamma- 
ray data (including early Fermi results) also put constraints on the dark matter decay 
explanation o f the PAMELA data, excluding some, but not all, of the relevant parame- 
ter space (e.g. Ishiwata et al. 201o[ Papucci &: Strumia 201ol ; 



Cirelli et al. 2010; Chen et al. 



20 id iMeade et al.lboiol : Izhang et al.lboiObl : IPohl & Eichlerlboioh . These studies typically 



utilize Fermi-LAT measurements of the diffuse Galactic and isotropic gamma-ray emission or 
Fermi sky maps and require that the expected decay signal from the Milky Way halo or the 
summed contribution from extragalactic dark matter not overproduce the observed gamma- 
ray fluxes. Limits on the dark matter lifetime from the o bservation o f dwar f galaxies with 
atmospheric Cherenkov telescopes have also been derived in lEssig et al.l (120091 ). Interestingly, 
these analyses indicate that current gamma-ray constraints on the decaying dark matter sce- 
nario are generically weaker than those on annihilating dark matter, with respect to the 
regions of parameter space providing an explanation to the cosmic-ray electron-positron 
anomalies. Moreover, as we review in the next section, several theoretical arguments imply 
that dark matter lifetimes in the range needed to explain the cosmic-ray ano malies are at least 



plaus ible, specifically in the context of supersymmetric grand unification (lArvanitaki et al. 
2009h . 



Beyond the interest in scenarios where decaying dark matter might be related to the 
cosmic-ray anomalies, several particle physics frameworks feature dark matter candidates 
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with a long, but finite, lifetime. We provide an overview of some of these models in the next 
section. The present study targets this general class of theories as well, and we consider 
here decay final states that can be regarded as representative for a wide range of scenar- 
ios. Searches for dark matter decay with gamma ray observations are therefore not only 
phenomenologically motivated, as they can test the dark matter interpretation of the ob- 
served cosmic-ray anomalies, but they also probe theoretically compelling extensions to the 
Standard Model of particle physics. 

While already significant, much can be done to improve the current limits on decaying 
dark matter models. In this work, we show that novel and significant constraints can be 
placed on the dark matter decay lifetime for both quark-antiquark decay final states as 
well as for leptonic final states giving a good fit the cosmic-ray anomalies using Fermi- 
LAT observations of isolated dark matter dominated objects like nearby cluste r s of g alaxies, 



dwarf spheroidal galaxies, and M31. In fact, the simulations of ICuesta et al.l ( 120101 ) which 
compared the gamma-ray flux from dark matter annihilation and decay in nearby large-scale 
structures as reconstructed from constrained cosmological simulations predict that nearby 
clusters and filaments are good targets for searches for and constraints on a signal from dark 
matter decay, and here we derive such constraints with recent Fermi observations. In §2, we 
discuss and introduce selected dark matter decay models, in §3 we discuss the selection of 
targets and the modeling of their dark matter content, and in §4 we present our results on 
the constraints on the decay lifetime. 



2. MODELS OF DECAYING DARK MATTER 

While stable on the scale of cosmological structure formation, the particle making up 
the universal dark matter might well have a finite lifetime, r. If this is the case, and if the 
dark matter particle decays into (visible) Standard Models particles, the decay process can 
be indirectly detected. Depending on the particle lifetime, non-thermal secondary particles 
resulting from dark matter decay can affect key "episodes" in the early universe such as the 
synthesis of light elements or the decoupling of the cosmic microwave background (CMB). 
In turn, light elemental abundances and the CMB spectrum can be used to set constraint s 



on particle models where the dark matter is unstable (see e.g. IChen fc Kamionkowskil 12004 ). 



There is no shortage of particle models where the universal dark matter is unstable, 
including unstable supersymmetric particles, moduli as dark matter particles, super- heavy 



dark matter candidates, axinos and sterile neutrinos ( IChen &: Kamionkowskil l2004j ). In- 
terestingly, particle candidates that would otherwise be completely un-detectable, such as 
gravitinos in supersymmetric extensions of the Standard Model, can in principle provide in- 
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direct clues on their nature if they decay ( IBerezinsky et al.lll99ll ). Among the decay particle 
debris, gamma rays and antimatter (such as energetic cosmic-ray positrons, anti-protons and 
anti-deuterons) are ideal signals for indirect dark matter detection. As such, gamma-ray and 
cosmic-ray antimatter data have been extens ively used to s et constraints on the lifetime of 
meta-stable particle dark matter candidates ( lBertonell201Cll ). 



The recent detection of an exceptionally large flux of high-energy cosmic-ray positrons 
(E e + > 10 GeV), compared to conventional models of Ga lactic secondary po sitron pro- 
duction, reported by the PAMELA space-borne experiment (jAdriani et al.ll2009l ) prompted 
great interest in possible non-standard sources. Another recently launched satellite experi- 
ment, Fermi-LAT , also showed the existence of a gentle but clear feature in the spectrum of 
cosmic-ray electrons-plus-positrons at energies of several hundreds of GeV, especially when 
suppl emented with eve n high er energy data from the H.E.S.S. atmospheric Cherenkov Tele- 



scope 



Aharonian et al.l ( 120081 ). While possibly pointing to nearby astrophysica l sources such 



as mature (characteristic age > 10 5 yr ) pulsars (see e.g. iHooper et al.l l2009t lYuksel et al. 



20091 ; |Profumoll2008l ; iGrasso et al.ll2009l ). the feature detected by Fermi-LAT and the rise in 
the PA MELA positron frac tion might share a common, more exotic origin: Galactic dark 
matter (jAdriani et al.l 120091 ) . A large number of studies have, in fact, demonstrated that 
particle models can be constructed to accommodate both the positron rise and the Fermi 
"feature" (for short, we will hereafter dub these two as the "cosmic-ray a nomalies" ) , invoking 
eithe r the pair-annihilation or the decay of dark matter particles (see e.g. iPapucci &: Strumia 
201dh . 



A dark matter annihilation interpretation of the cosmic-ray anomalies requires a very 
large pair-annihilation cross section, one or more orders of magnitude larger (depending 
on the dark matter particle mass) than the one that would ex plain via thermal decoupling 
the observed universal matter density ( iBergstrom et al.ll2009l ). The attractiveness of the 
thermal relic dark matter scenario must therefore be traded-off, within this scenario, for 
large annihilation rates in today's cold universe, unle ss a very strong velocity d e pendence 



in the annihilation cross-sect ion is invoked, (see e.g. lArkani-Hamed et al.l l2009t llbe et al. 



20091 ; ICyr-Racine et al.l 120091 ) . Dark matter decay as the source of the cosmic-ray lepton 



anomalies, on the other hand, requires the fine-tuning of the decay time-scale r, but is 
otherwise a phenomenologically perfectly viable scenario. The required lifetimes are on the 
order of a few 10 26 s, for particle masses in the TeV range. 

In the broader picture of particle model-building, several scenarios predict that the 
symmetry that preserves the stability of the dark matter particle on scales much larger than 
the age of the universe is violated by sufficiently suppressed higher-dimension operators (i.e. 
operators with mass-dimension greater than 4) that mediate the subsequent, slow particle 
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decay. Similarly to the generic prediction of proton decay, in the context of supersymmetric 
grand unification, operators with mass dimension 6 are expected to make supersymmetric 
dark matter unstable, with a natural time-scale 



87rMG UT /m 5 



DM; 



(1) 



where Mqtjt is the energy scale of grand unification interac t ions, as recently pointed out by 
Nardi et al.l (120091 ) and, subsequently, by lArvanitaki et al.l ( 120091 ) . Interestingly this yields 
the prediction 



T 



10 



27 



sec 



/ TeV V 

\ m DM / 



Mi 



GUT 



m DM / \2xl0 16 GeV, 
which is precisely in the ballpark needed to explain the cosmic-ray anomalies. Notice that 
GUT-suppressed dimen sion-5 operators would produce lifetimes on the order of a second 
(lArvanitaki et al.l 120091 ). and must thus be suppressed to avoid dark matter decaying long 
before structure formation. 

In addition to grand unification, another compelling reason to expect dark matter to 
decay with a lifetime on the order of 10 26 — 10 27 s is provided by models, such as those eluci- 



dated in the analysis of llbarra fc Tranl (120091 ). where the lightest supersymmetric particle is 
a gra vitino and .R-parity is (very weakly) broken (see also iBertone et al.ll2007l ; llbarra &: Tran 
20081 ). In this class of models, not only is the supersymmetric gravitino problem naturally 
solved, but the correct baryon asymmetry can naturally be generated by i?-parity violating 
interactions via leptogenesis. Another scenario that might in principle combine the origin 
of both the baryonic matter and the dark matter via leptogenesis, and provide a mecha- 
nism to_expl_a2n_rj^tnno masses and mixing is that of sterile neutrinos (for a review see 



e.g. iBoyarsky et al.ll2009l ). Recent analyses also pointed to hidden-sector gauge bosons and 
gauginos as natural particle frameworks where TeV-scale particle dark matter decays with 
lifetimes on the order of 10 26 s. A scenario featuring leptonic decay final states and a suffi- 
ciently long-lived dark matter as a consequence of a global symmetry broken by instanton 
op erators origin a ting f rom a non-Abelian gauge group in the dark sector was also discussed 



m 



Carone et al. (boioh 



If the decay of particle dark matter is to explain the cosmic-ray lepton anomalies, the 
re quired dominant decay fina l state must be close to one (or more) muon pairs, as shown e.g. 



in 



Papucci fc Strumial (j2010l ). Nevertheless, among many of the scenarios mentioned above, 



the decay products will likely follow hadronization chains that are typically well reproduced 
in their major features by de cay final states such a s a quark-antiquark pair (even though 
for instance in the scenario of llbarra fc Tranl (120091 ) the final state might be that resulting 
from W or Z decay, the hadronic decay modes would still pro duce a gamma-ray spe ctrum 
qualitatively similar to that of a quark-antiquark pair; see also ICiafaloni et al.l (120101 ) for a 
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discussion of weak corrections to annihilation or decay final states). For this reason, in the 
present study we will focus on two specific decay final states: fi + fi~ as representative of those 
scenarios that offer a particle dark matter decay interpretation to the cosmic-ray anomalies; 
and bb for more conventional decaying dark matter scenarios. 

For these two decaying dark matter final states, we will consider in detail the pro- 
duced gamma-ray spectrum as a functi on of the dark matter par ticle mass, and model 
this spectrum with the DMFIT package (IJeltema &: Profumoll2008al ) as implemented in the 
Fermi-LAT data analysis pipelin^j]. The consideration of different decay final states will 
be coupled to the choice of the best targets, as one particular target does not necessarily 
gi ve the strongest con s traint s for all dark matter models . An important ef f ect, as shown 



m iColafrancesco et al.l (120071 ) , IJeltema &: Profumol (l2008bl ) , and lAbdo et al.l ( 1201 Oq ) , is the 
diffusion of the high-energy electrons and positrons produced by dark matter decay. For 
sufficiently massive dark matter candidates, in fact, these particles copiously up-scatter low- 
energy back ground photons (inv erse Compton - IC - emission) into gamma rays detectable 
by the LAT ( Zhang et al. 2010bl Jah . The IC photon flux is however strongly dependent upon 
the electron diffusion and leakage out of the observed target, an effect that impacts both 
the spectrum and the predicted flux-intensity. We discuss this and other target- dependent 
considerations in the next section. 



3. TARGETS FOR DARK MATTER DECAY SEARCHES 

A number of astrophysical targets are potentially good sources to search for a signal from 
dark matter annihilation or decay, including the Galactic center, the Milky Way halo, nearby 
galaxies and dwarf spheroidal galaxies, nearby clusters of galaxies, and the summed contri- 
bution from unresolv ed extragalactic d ark matter (for a review in the context of the LAT 



performance, see e.g. iBaltz et al.l 120081 ) . Underpinning which is the "best" target depends, 
however, on the type of signal (annihilation or decay), on our understanding of the dark 
matter density in different objects, on the gamma-ray backgrounds, and on the dominant 
gamma-ray emission mechanism for a given dark matter particle model (for instance prompt 
gamma-ray production versus secondary IC emission). Extragalactic objects like clusters of 
galaxies and Local Group dwarf galaxies are particularly attractive targets for dark matter 
searches as they are highly dark matter dominated, isolated objects, and many of them lie 
in regions of low diffuse gamma-ray background, being located at high Galactic latitudes. In 
addition, with the exception of M31, the objects we discuss do not themselves host detectable 



1 http://fermi.gsfc.nasa.gov/ssc/ 



-7- 



gamma-ray emission, making them clean targets of study to constrain emission from dark 
matter. 

Previous studies looking at constraints on dark matter models usi ng Fermi observa- 



tions of dwarfs and clusters have focused on dark matter annihilation (Ackermann et al. 



2010bl ; lAbdo et al.ll2010bl ). while constraints on dark matter decay with Fermi have primar- 
ily compared the observed diffuse gamma-ray backgrounds to the total summed emission 
expected fro m dark matter decay in the Milky Way halo or from unresolved extragalactic 



dark matter (llshiwata et al. 



2010: Meade et al. 



2010; 



Papucci fc Strumiall^uiui; imreii 
boiObi IPohl fc Eichlerlboioh 



20101 ; IZhang et al. 



2010; 



Cirelli et al. 2010; Chen et al. 



In this paper, we con- 
strain dark matter decay models using, instead, isolated extra-galactic dark matter sources, 
specifically nearby clusters, dwarf spheroidal satellites of the Milky Way, and the nearby 
galaxy M31 (the Andromeda galaxy). In some cases, these targets provide significantly 
stronger constraints on the dark matter decay lifetime than have been previously published , 
particularly for the case of nearby clusters (see also the simulations of ICuesta et al.ll2010l ). 
which we discuss in the next section. 



3.1. CLUSTERS OF GALAXIES 



Clusters of galaxies are the most massive collapsed objects in the present-day Universe 
and are highly dark matter dominated. Apart from the central bright active galactic nuclei in 
a couple of clusters, clusters have thus far not been detected in ga mma-rays, a fact that allow s 
one to place limits on t heir possible dark m atter related emission ( Ackermann et al. 2010bl 



Aharonian et al. 2009; Aleksic et al. 



2010 ). We consider here the same sample of six, nearby 



clusters of galaxies as we studied in I Ackermann et al.l (j2010bl ) . These clusters are among the 
best candidates for dark matter searches and have published upper limits on their potential 
gamma-ray emission from dark matter (annihilation) based on the first 11 months of Fermi- 
LAT observations for two illustrative (annihilation) final states, bb and ^ + fi~, and for a grid 
of dark matter particle masses. Remarkably, the expected gamma-ray spectral shape from 
dark matter decay for a given particle mass is the same as the gamma-ray spectrum from 
annihilation for a particle of half the mass with the same assumed final state (although certain 
fina l states can only be available to decay pr ocesses, such as e.g. W plus charged leptons, see 
e.g. iPalomares-Ruiz fc Siegal-Gaskinsll2010r) . and we use the corresponding flux upper limits 
from lAckermann et al.l (l2010bt ) to take into account the spectral dependence of the broad- 
band Fermi-LAT upper limits. Clusters are particularly good targets when considering decay 
to leptons, as is the case for models fitting the cosmic-ray anomalies, because in this case 
much of the gamma-ray emission stems from IC scattering of background radiation by the 
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energetic e + e~ produced (lAckermann et al.ll2010bl ; IColafrancesco et al.ll2006l ). In clusters, 
e + e~ lose energy through IC much f aster than they can diffuse out of the system, unlike 
in sm aller systems like dwarf galaxies (jJeltema fc Profumoll2008bl ; IColafrancesco et al.ll2007l . 
2006h . 



The Fermi-LAT effective area is a strong function of energy and increases by nearly 
an order of magnitude from 100 MeV to 1 GeV, increasing more slowly at energies above 1 
In addition, the LAT PSF decreases with increasing energy (effectively reducing the 
background for detection) and the gamma-ray backgrounds are relatively soft. Given the 
expected spectral shape for the dark matter decay models we consider, these facts mean that 
the tightest constraints come from the Fermi data around ~ 1 — 2 GeV (for more detail see 



e.g. lAckermann et al.l l2010bU afl I. At these energies the Fermi-LAT Point-Spread Function 
(PSF) is around on e degree, a radius which co ntains most of the virial mass of the clusters we 
consider. Following lAckermann et al.l ( j2010bl ). we assume a n NFW profile for the cluster dark 
matter density distribution with cl uster masses taken f rom lReiprich fc Bohringerl (120021 ) and 
the concentration-mass relation of iBuote et al.l (120071 ). The normalization of the expected 
gamma-ray flux from dark matter decay Jd is 



J d = [ dSl [ pnu(l)dl(ip), (2) 
J An ji.o.s. 

where we integrate the dark matter density pdm along the line of sight (l.o.s.) in the direction 
t/> over a solid angle AQ corresponding to a radius of one degree. The Jd values for our cluster 
sample are listed in Table 1 along with the masses, radii, and concentrations we employ. The 
choice of density profile and integration radius only matters at the level of a factor of a few 
in the derived Jd since, as already mentioned, the size of the Fermi PSF includes much of 
the cluster virial mass (see column 6 of Table 1) and since compared to annihilation, which 
is proportional to density squared, decay is relatively insensitive to the inner dark matter 
density profile. 



3.2. DWARF SPHEROIDAL GALAXIES 



Local Group dwarf spheroidal galaxies are v ery nearby objec ts which, in some cases, are 
even more dark matter dominated than clusters (IWolf et al.ll2010l ). Similar to clusters, dwarf 



2 http://www-glast.slac.stanford.edu/software/IS/glastJat_performance.htm 

see also the LAT point source sensitivity for an assumed power law spectrum here: http://www- 
glast.slac.stanford.edu/software/IS/glastJat_performance.htm 
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spheroi dal galaxies have not been detected in gamma-rays ( lAbdo et al.ll2010bl ). lAbdo et al. 
( I2010bl ) places limits on dark matter annihilation models based on the non-detection of 
eight dwarf spheroidal galaxies in the first 11 months of Fermi-LAT data; we will consider 
this sample of eight dwarfs in terms of limits on dark matter decay employing again the 
appropriate published Fermi-LAT flux limits for assumed bb and fi + fi~ final states and a 
grid of particle masses. 

For dwarf spheroidal galaxies, the dark matter halo mass can be estimated from the 
line-of-sight stellar velocity dispersion. Dwarf masses are most strongly constrained near 
the half light radius as the mass at the half light ra dius, M1/2, is re latively insensitive 



to the unknown stellar v elocity dispersion anisotropy ( IWolf et al.ll2010l ) (see however also 
Amorisco fc Evansll2010l ). Outside of this radius the mass and dark matter density profile 
are less wel l constrained (though the use of cold dark matter (CDM) priors reduces this 
uncertainty iMartinez et al.ll2009l ). The half light radii of the dwarf spheroidal galaxies we 
consider are smaller than the Fermi-LAT PSF at our energ ies of interest (se e column 5 of 
Table 2), so we will conservatively consider M1/2, taken from [Wolf et al.l ( 120101 ). as the lower 
limit on the dark matter mass probed by the Fermi observations. We list the corresponding 
values of as well as M1/2 and ri/2 in Table 2. For comparison, in Figure H] we show the 
constraints we would get f or a range o f plau sible dwarf virial masses assuming the density 
profile given in Figure 3 of IWolf et al.l ( 120101 ) and taking the corresponding mass within one 
degree. 



3.3. M 31 



As the nearest large galaxy, M31 (the Andromeda galaxy) is also an excel lent target for 



dark m atter decay searches, as envisioned e.g. for the case of sterile neutrinos in lWatson et al. 



( 120061 ) . A significant detection of M31 has recently been reported in two yea rs of Fermi-LAT 
data (IBechtol fc the Fermi-LAT Collaboration! |2010| ; lOgelman et al.l |2010| ). The gamma- 
ray luminosity of M31 is consistent with predictions for the luminosity from cosmic ray 
collisions with the ISM a nd is likely not due to dark matter. We will conservatively take 
the detected flux plus 2a (IBechtol fc the Fermi-LAT Collaborationll2010l ) as the upper limit 
on the possible dark matter induced gamma-ray emission from M31 in order to place limits 
on the dark matter decay lifetime, though in all likelihood any dark matter signal from 
M31 is significantly fainter. We consider both the flux resulting from the assumption that 
the gamma-ray emission is point-like and from the assumption of an extended source with 
emission following the M31 infra- red surface br i ghtne ss. To find Jj, we employ the NFW 
dark matter density profile from iGeehan et al.l ( 120061 ) based on fits to surface brightness, 



velocity dispersion, and rota tion curve data (see Table 1). We also consider the isothermal 
profile of iKerins et al.l ( 1200 ll ). for comparison. 



4. CONSTRAINTS ON DECAYING DARK MATTER MODELS 

In this section we present our results on the limits on the decay lifetime for dark matter 
particles as a function of mass from gamma-ray observations of clusters ( §4.11) . of dwarf 
spheroidal galaxies ( §4.2p and of M31 ( §4.3[) . As a rule of thumb, it is simple to give an 
estimate for the expected Fermi-LAT limits for the case of a standard decaying WIMP: 

1Q- 10 cm-^ec' 1 ^ / J d \ / 100 GeV \ / Ky\ 

^ J Vl0 18 GeVcm-V \ m DM J \~10 ) ( ^ 

where 7 is the limit on the gamma-ray flux from the observed target, and iV 7 the integrated 
number of photons per decay within the relevant Fermi-LAT energy range (this number 
depends on the decay final state as well as on the decaying particle mass m DM ). 



r < 8 x 10 25 



sec 



4.1. CLUSTERS OF GALAXIES 



Figure 1 shows the limits we derive for the possible dark matter decay lifetime assuming 
a bb final state (left panel) and final state (right panel) derived from the 95% confi- 

de nce level upper limits on the gamma-ray flux from nearby clusters of galaxies presented 



m 



Ackermann et al.l (l2010bl ) for 11 months of Fermi-LAT observations. For comparison, we 



show previous constraints on the decay lifetime derived from Fermi all-sky maps and consid- 
ering the dark matte r contribution from b oth the M ilky Way halo and unresolved extragalac- 



tic dark matter from IZhang et al. fhoiobl) (see also llshiwata et al.ll2010t iPapucci fc Strumial 



2010: Cirelli et al. 2010; Chen et al. 2010; Meade et al 



2010 



Pohl fc Eichleril2010l ) 



The shape of the constraints on the lifetime versus mass can be understood by consid- 
ering the combined effect of how the gamma-ray flux depends on the decaying particle mass 
and of how the flux sensitivity depends on the resulting gamma-ray spectrum. Increasing 
the particle mass mDM, the particle number density, and hence the number of decays, de- 
creases as wdm oc 1/mDM- However, the number of photons per decay increases with mDM - 
more energy density is dumped into the decay products overall as the decaying particle mass 
is increased. The photons relevant to derive constraints from Fermi-LAT data, however, 
are restricted to a specific energy range. Particularly when secondary photons from Inverse 
Compton scattering are included, deriving a simple functional form for the integrated num- 
ber of photons A" 7 within the Fermi-LAT energy range from the decay of a particle of mass 
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7«dm versus the particle mass is non-trivial. Suffice it to say that, neglecting IC photons iV 7 
grows with ttj-dm slower than linearly. This implies that, when using a fixed flux upper limit 
(as we do for the case of M31, see below) we expect the limits on the dark matter lifetime 
to decrease as m^^, with < a < 1. 

A more subtle effect occurs if the flux limits are calculated themselves as a function of 



mass, a s is done in the dedicat ed Fermi-L AT data analyses presented in lAckermann et al. 



( 12 01 Obi ) and I Ab do et al.l ( 12 01 Obi ) where the simulated gamma-ray spectra from annihilation of 
particles into given final states were employed. The flux upper limits improve (i.e. decrease) 
with the particle mass, since most of the photons are produced at increasingly large energies 
where the instrum ental point-spread function and effective area improve (see e.g. Fig. 2 in 
Abdo et al.ll2010bl ). This effect is even more acute when photons from IC are accounted for 
- in this case for large enough masses m DM > 100 GeV, photons start to enter the Fermi 



LAT e nergy range and to dominate the total photon count (see Fig. 2 in lAckermann et al. 



2010b|). As a result of the combination of the flux limits improving with tbdm and of the 
photons counts decreasing with mDM, we expect a maximum, at some intermediate value of 
the particle mass, for the constraints on the lifetime. This indeed what we find, unless, as 
in the case of M31, we simply employ one fixed flux limit value. 

In the bb panel, the yellow band to the left shows a parameter spac e range corresponding 
to the oretical predicti o ns ba sed on the gravitino dark matter model of iBertone et al.l (120071 ) 
and of llbarra fc Tranl ( 120081 ). where we posit that the gravitino lifetime r 3 / 2 goes as 



T3/2 



3.8 x 10 27 s 



10 



-8 



m 



3/2 \ 



-3 



10 GeV 7 



with m 3 / 2 the gravitino mass, and the mo del- dependent parameter U i 



grav 



10" 



(4) 



fin the 



figure, the left boundary corresponds to U gTSCV = 2 x 10~ 8 and the right boundary to U gTajv = 
0.5 x 10~ 8 ). The right-most yellow band, in the bb panel, indicates a range for the mass-scale 
Mnp of a "New Physics" dimension-6 operator that would mediate dark matter decay in the 
context of supersymmetric grand unification, according to the rate outlined in Eq. ([1]) with 
M GUT = M NP . We take as an indicative range 10 16 < M NP /GeV < 4 x 10 16 . 

In the right panel, we show decay lifetimes giving a good fit to the PAMELA positron 
excess and the combination of the PAMELA position fraction and the feature observed in 
the e + + e~ spect rum measured by the F ermi-LAT (for details on the determination of 
these contours, see Papucci fc Stmmia] 2010 ). We also shade in yellow the more "optimistic" 
limits calculated by IZhang et al.l (l2010bl ). where a specific, albeit rather arbitrary, Galactic 
gamma-ray diffuse model is assumed and subtracted off of the Fermi data. 



It can be seen from Figure 1 that clusters of galaxies give significantly stronger con- 
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straints on the dark matter decay lifetime, both for a bb and a final state and for a 

large range of particle masses, than have previously been derived from gamma-ray observa- 
tions. Specifically, for a bb final state the decay lifetime is limited to be above ~ 5 x 10 26 
sees for particle masses in the very broad mass range between 10 GeV and 1 TeV, while 
previous limits were about a factor of five worse worse. Limits on specific particle theories 
with meta-stable dark matter, such as those shown in the figure with shaded yellow contours, 
also improve accordingly. 

For a final state (right panel), the non-detection of clusters by Fermi excludes 

essentially the entire best fit region for dark matter decay models fitting the combination 
of the PAMELA positron fraction and the Fermi-LAT e + + e~ spectrum and much of the 
region fitting the P AMELA data alone. Notice that this is in line with the predictions of 



Pinzke et al.l (120090 . where the effect of substructures was also extensively discussed. Only 
models where the dark matter particle is significantly lighter than 2 TeV seem to pass the 
test of gamma ray emission from galaxy clusters. These results, therefore, cast doubt on the 
dark matter decay explanation of the rising positron fraction at high energies, particularly if 
this is taken jointly as an explanation of the Fermi-LAT electron-positron spectrum. Similar 
to dark matter annihilation, the strongest constraints come from gamma-ray observations of 
nearby groups or poor clusters like Fornax, M49, and NGC4636, but all of the clusters in our 
sample lead to an improvement in the dark matter decay constraints. We note from Table 
1 that the ranking in terms of the brightest signal from dark matter of the various clusters 
is quite different for the cases of decay and annihilation, although the brightest object is in 
both cases the Fornax cluster. Notice that some uncertainty is present in the determination 
of the coefficient Jd, as we discuss below. 

In Figure [2j we compare the constraints on the dark matter decay lifetime determined 
above (blue lines) using an X-ray derived dark matter density profile to alternate choices for 
the density profile or cluster mass as an indication of the potential systematic error in our 
constraints for the Fornax cluster (left panel), which has the highest Jd of the cluster sample, 
and for the Coma cluster (right panel), which has the highest mass of the cluster sample. 
We only explicitly show decay to but the same relative offsets hold for other final 

states. For the Forna x cluster, we plot t he lim its for two alternate estimates of the cluster 



mass: the mass from iDrinkwater et al.l (120011 ) based on the velocity dispersion of cluster 



member galaxies and the mass from the independent X-ray analysis of iPaolillo et al.l (120021 ) 
based on ROSAT HRI observations (mass within one degree given from their published mass 
profiles). The dark matter density profile in the central Fornax galaxy, NGC1399, has also 
been investigated based on the radial velocities of globular clusters (GC), but GC dynamics 
(particularly the red GC) probe primaril y the central galaxy its elf and do not reliably give 



the potential of the cluster as a whole (jSchuberth et al.l |201Q| ). The mass of the Coma 
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cluster, the most massive and most distant cluster in our sample, has also been estimated 
based on weak lensing of background galaxies. In the right panel of Figure |2l we show the 
change in the constraints on the dark matter decay lifetime for three different weak lensing 



estim ates of the Coma cluster mass (IKubo et al.l 120071 ; iGavazzi et al.l l2009t lOkabe et al. 



20101 ) as well a s for the cluster mass estimated based on the velocity dispersion of the 



cluster galaxies (ILokas fc Mamonl 120031 ) . For both Fornax and Coma, the constraints for a 
variety of mass estimates change by less than a factor of two from the nominal constraints 
shown in Figure 1. Even the most conservative Fornax constraints represent a significant 
improvement over previous limits and exclude large portions of the parameter space that 
could explain the cosmic ray anomalies. Fo r our nominal constraints we employ the mass- 
concentration relation of iBuote et al.l (120071 ). because this relation was derived based on X- 
ray cluster observations; however, we not e that for a variety of mass-concentration relations 



from both simulations a nd observations ( jZhao et al.l 120091 ; iMaccio et al.l 120081 ; iKlypin et al. 



2010l ; lEttori et al.l 120101 ) the Jd of Fornax would change by less than 20%. 



A full accounting for systematic effects for all clusters in the sample is not statistically 
not feasible (the observations not being directly comparable), nor do we wish to debate the 
merits of one set of observations over another. However, Figure [2] leads us to conclude that 
our constraints are robust within roughly a factor of two, and that they are not systematically 
biased high or low. As the Fermi non-detection of most of the clusters in the sample excludes 
a significant fraction of the decay lifetime and particle mass combinations needed for a joint 
dark matter decay explanation to the cosmic-ray anomalies, our results lead to the conclusion 
that this possibility appears disfavored. 



4.2. DWARF SPHEROIDAL GALAXIES 



Figure [3] shows the limits on the dark matter decay lifetime for nearby d warf spheroidal 
galaxies using the Fermi gamma-ray limits published in lAbdo et al.l (j2010bl ). For the dwarf 
constraints with a yU + /i~ final state (right panel), we conservatively only include gamma- 
ray emission from final state radiation and do not include I C emission, given the un 



known effects of diffusion on e + and e in dwarf scale systems (jJeltema fc Profumoll2008b 



Colafrancesco et al.l 120071 ) . Below we discuss how these constraints would improve with the 
inclusion of IC gamma-ray emission and specific assumptions for cosmic-ray diffusion. Com- 
pared to clusters, the non-detection of dwarf spheroidal galaxies in gamma-rays gives weaker 
constraints at least for the conservative limits we derive assuming a mass of M\/2 and no 
IC emission. However, at low masses the constraints still improve upon previous limits from 
the literature, in particular those we derive from the Ursa Minor dwarf galaxy. 
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If we relax our conservative assumptions on gamma-ray emission from high-energy muon 
pairs, the constraints shown in the right panel of Figure [3] become considerably stronger, 
as shown in Figure H] for Ursa Minor (we employ this specific dwarf galaxy as it gives the 
strongest constraints among the dwarf sample). In the left panel of Figure IU we consider the 
limits on the dark matter decay lifetime if we include the expected IC gamma- ray emission 



for a u + U decay final state for different assumed diffusion coefficients. Following I Ab do et al. 



( 12010bl ). the blue lines show the constraints for diffusion coefficients of -Do = 10 28 cm 2 s _1 
(dot-dashed line) and Dq = 10 29 cm 2 s _1 (dashed line) which bound the values found for the 
Milky Way. In both cases, the energy dependence of diffusion is assumed to have a power-law 

1 /3 

dependence given by D(E) = D ( x Gev ) • ^ ne rec ^ nne snows the constraints in the limit 
of no diffusion (i.e. in the limit where the electrons and positrons loose all of their energy 
via IC prior to diffusing away from the system). The inclusion of IC emission significantly 
strengthens the constraints from Ursa Minor, but for reasonable diffusion coefficients these 
are still weaker than previous limits on the decay lifetime. 

In addition for our nominal constraints in Figure |3j we conservatively assumed dwarf 
masses equal to M\/2 even though the dwarf half-light radii are smaller than 1 degree. In the 
right panel of Figure HI we compare the constraints for Ursa Minor for different assumptions 
on the mass profile. For illustration we show the limits on the decay lifetime for a 
final state with a diffusion coefficient of D = 10 28 cm 2 s -1 , but of course the same relative 
offset also holds for other cases. Employing the best-fit NFW profi l e para meters from fits 
to the stellar velocity data for Ursa Minor published in lAbdo et al.l ( l2010bl ) to calculate Jd 
for a radius of 1 degree yields the d o t-das hed line in Figure [4] (right panel), a factor of a few 
stronger than for M\i2- IWolf et al.l ( 120101 ) find that all of the Milky Way dwarf spheroidals 
are consistent with having the same total halo mass (within the scatter) of around 3 x 10 9 
M for a ACDM cosmology. Looking for example at their Figure 3, for a range of two orders 
of magnitude in total halo mass between 3 x 10 8 M Q and 3 x 10 10 M Q and for an NFW profile 
consistent with a WMAP5 cosmology, would give a mass within 1 degree for Ursa Minor 
between a few times 10 7 M Q and a few times 10 8 M . The decay constraints for this range 
of masses is reflected in the two dashed lines in Figure HI where we show the results for the 
specific values of M(r < 1°) = 3 x 10 7 M and 3 x 10 8 M Q . Depending on the assumptions 
on mass and diffusion, the Ursa Minor constraints become comparable to or stronger than 
previous limits. 
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4.3. M 31 



Although faintly, the nearby spiral galaxy M31 is detected at gamma-ray frequencies, as 
recen tly reported by the Fermi-LAT Collaboration (IBechtol &: the Fermi-LAT Collaboration 
2010l ). The detected gamma-ray emission is compatible with a point-source, although a 
slightly higher test-statistics w as reported if a spatial template based on the infra-red sur- 
face brightness was employed (IBechtol Sz the Fermi-LAT Collaboration! 120101 ). Here, to set 
constraints on the emission from dark matter decay we simply consider the 2-cr upper limit 
on the observed emission, conservatively summing the statistical and systematic uncertain- 
ties quoted by the Fermi-LAT Collaboration. A more in-depth study would resort to either 
attempting to subtract the astrophysical emission from M31 as expected from cosmic-ray 
interaction with the inter-stellar medium (in analogy to the diffuse gamma-ray emission from 
our own Milky Way), or to utilize, for instance, a region of interest that cuts out the gamma- 
ray emitting region (e.g. an annulus around the M31 central region). We defer such a study 
for a forthcoming analysis. 

M31 is an extremely interesting object to search for a signal from dark matter decay, 
as was for instance suggested, for the case of sterile neutrinos (producing a monochromatic 
line at X-ray frequen cies from the two-body decay into an active neutrino and a photon) by 
Watson et al.l ( 120061 ). To estimate the object's brightness in gamma r ays from dark matte r 
decay, we employ the NFW dark matter density profile proposed by iGeehan et al.l ( 120061 ): 
Table 1 shows that we obtain a Jd value (1 degree radius) for M31 which exceeds by almost 
a factor 2 the largest normalization factor for galaxy clusters, and by more than an order of 
magnitude the largest one we find among dwarf galaxies. 

Despite being a potentially extremely bright source of gamma rays from dark matter 
decay, the constraints we derive from M31 are not as strong as those we found for clusters 
- the reason evidently being that M31 is actually detected as a gamma-ray source, and we 
conservatively do not attempt to subtract off any astrophysical background. The left panel 
of Figure illustrates the constraints we obtain for a bb decay final state. The black line 
refers to the NFW profile with the gamma-ray flux obtained with the assumption of a point- 
lik e source and Jd as l isted in Table 1. The green dashed line assumes the isothermal profile 
of iKerins et al.l ( 1200 ll ) and the same gamma-ray flux limit - the closeness of the two lines 
illustrates the small uncertainty that stems from using two entirely different functional forms 
for the dark matter density profile. Finally, for the blue dotted line we employ the extended 
source Fermi-LAT flux and assume that the relevant mass coincides with the virial mass of 



M31 (which we set to M 2 oo — 6.8 x 1O 11 M , from lGeehan et al.ll2006[ ). The limits we obtain, 
while an improvement over what we get from dwarf galaxies and over previously available 
limits, are a factor of a few worse than those implied by the non-detection of galaxy clusters 
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for the bb final state. 

The right panel of Figure shows our results for the \x + \x~ final state. We compare again 
with the regions favored by a dark matter decay interpretation of the cosmic-ray anomalies 
(green and red shaded regions) and with the current constraints (grey region). The black 
lines refer to the NFW plus point-like source assumption, while the red lines refer to the 
case where we employ the virial mass and the extended source gamma-ray flux. Given the 
relevance of diffusion for cosmic-ray electrons and positrons in the Milky Way, we do expect 
leakage of secondary e + e~ out of M31 as well, and we model this process again in analogy 
to our own Galaxy. We show with solid lines our results for the relatively small diffusion 
coefficient Do = 10 28 cm 2 /s, and with dashed lines the more conservative case of Do = 10 29 
cm 2 /s (as for dwarf galaxies, we set the rigidity dependence power to 1/3). In either case, 
the limits we obtain do not improve over previous constraints nor do they approach the 
phenomenologically relevant regions for the cosmic-ray anomalies. 

5. DISCUSSION AND CONCLUSIONS 

In this paper, we analyzed how early Fermi limits on the gamma-ray flux from clusters 
of galaxies, from nearby dwarf galaxies and the Fermi detection of M31 impact the possible 
lifetime of a meta-stable particle dark matter. We showed that limits from dwarf galax- 
ies and from M31 improve over previous constraints, for certain ranges of masses and final 
states. The best limits are however derived in the case of galaxy clusters, where we obtain 
an improvement of about a factor of five on the allowed decay lifetime for hadronic decay 
final states for a wide range of particle dark matter masses between 10 GeV and 1 TeV. 
Furthermore, models where dark matter decays into muon pairs and explains the so-called 
cosmic-ray anomalies (the positron excess reported by PAMELA and the feature observed 
by Fermi in the electron-positron spectrum) are found to over-produce gamma rays from 
clusters of galaxies compared to Fermi observations for large regions of the best-fitting pa- 
rameter space. While we cannot completely exclude the possibility, our results imply that 
the dark matter decay interpretation of the cosmic-ray anomalies is in tension with gamma 
ray observations, especially if both the Fermi and PAMELA anomalies are jointly explained, 
and this scenario may be ruled out completely if Fermi does not detect gamma-ray emission 
from clusters in the next few years. Fermi-LAT continues to take data in all-sky survey 
mode, and in the future longer exposures and a better understanding of the astrophysical 
backgrounds will continue to improve constraints on dark matter models and perhaps even 
reveal a signal from this elusive component of the universe. 
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Table 1. Data for Clusters and M31 



Name 


Distance 


M 200 


C200 


-R200 


Tldeg/ Rvir 


J d 


Ranking for 




Mpc 


1O 14 M 




Mpc 




10 18 GeV cm" 2 


DM Annihilation 


Fornax 


18.9 


1.42 


8.9 


1.35 


0.24 


18.4 


1 


Coma 


94.9 


19.38 


5.7 


3.22 


0.50 


16.6 


5 


Centaurus 


42.2 


3.78 


7.5 


1.87 


0.39 


13.7 


4 


M49 


18.1 


0.65 


10.2 


1.04 


0.30 


11.1 


2 


AWM7 


70.4 


6.08 


7.0 


2.19 


0.55 


10.2 


6 


NGC4636 


15.2 


0.35 


11.4 


0.85 


0.31 


8.88 


3 


M 31 


0.77 


0.0068 


22.0 


0.18 


0.075 


36.4 





Note. - For clusters, M 2 qo and -R200 are taken from iReiprich fc Bo hringer (2002) and the 



concentrations are calculated from the concentration-mass relation of 



BuoteetaL 



(120071 ): for 

M31 these quantities are taken from lGeehan et all (120061 ) . Column 6 gives the ratio of the radius 
which subtends a projected angle of 1 degree (i.e. the radius within which J d is calculated) to 
i?2oo- The Jd factors, listed in column 7, for the clusters under consideration and for M31 are 
calculated as described in §3. For comparison, column 8 lists the relative ranking of the same 
objects in terms of their normalizations for a signal from dark matter annihilation. 
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Table 2. Data for Dwarf Spheroidal Galaxies 



Name 


Distance 


Mi /2 


ri/2 




■h 


Ranking for 




kpc 


1O 7 M 


pc 




10 18 GeV cm" 2 


DM Annihilation 


Ursa Minor 


66 


5.56 


588 


2.0 


1.50 


2 


Ursa Major II 


30 


0.79 


184 


2.8 


1.03 


3 


Bootes I 


62 


2.36 


322 


3.4 


0.72 


5 


Sextans 


86 


3.49 


1019 


1.5 


0.55 


7 


Fornax 


138 


7.39 


944 


2.6 


0.45 


8 


Draco 


76 


2.11 


291 


4.6 


0.43 


1 


Sculptor 


79 


2.25 


375 


3.7 


0.42 


4 


Coma Berenices 


44 


0.197 


100 


7.7 


0.12 


5 



Note. — The distances to the dw arf spheroidal gala xies are taken from lAbdo et al.l (j2010bl ). 
while Mi/2 and 7*1/2 are taken from IWolf et al.l (120101 ). Column 5 gives the ratio of the radius 
which subtends a projected angle of 1 degree (i.e. the radius within which Jd is calculated) 
to 7*1/2. Column 6 lists the Jd factors for dwarf galaxies under consideration calculated as 
described in §3. For comparison, column 7 lists the relative ranking of the same objects in 
terms of their normalizations for a signal from dark matter annihilation. 
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Fig. 1. — Constraints on the decay lifetime as a function of mass from clusters of galaxies for 
a bb final state (left panel) and a final state (right panel) based on the 95% confidence 

level upper limits on the gamma-ray flux from 11 months of Fermi-LA T observations. The 
grey shaded region shows an example of previous constraints from IZhang et al.l (J2010bJ) 
which consider the expected diffuse dark matter decay signal from the Galactic halo and 
unresolved extragalactic dark matter; the yellow sh aded region in the right panel shows 
the more optimistic limits from IZhang et al.l (j2010bl ) obtained after subtraction of a model 
for the Galactic diffuse emission. In the left panel, the yellow bands show the regions for 
theoretically motivated models for gravitino decay (leftmost yellow band) and for theories 
with a dimension-6 operator at the GUT scale mediating dark matter decay (rightmost 
yellow band), see §4.1 for details. In the right panel, we also show the regions of parameter 
space fitting the observed cosmic-ray anomolies (green: PAMELA positron fraction only, red: 
the combination of the PAMELA position fraction and the feature o bserved in the e + + e~ 
spectrum measured by the Fermi-LAT, see lPapucci fc Strumiall2010l ). 
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Fig. 2. — Constraints on dark matter decay lifetime as a function of particle mass for 
alternate determinations of the mass of the Fornax (left panel) and Coma (right panel) 
clusters. The constraints are shown for decay to but the same relative offsets hold for 

other final states. The nominal constraints based on the X-ray determined masses within 1 
degree are reproduced from Figure 1 by the solid blue lines (see §3.1). Left: The red lines 
show the constraints if we employ, instead, the dark matter mass of Fornax (radius of 1 



Drinkwater et al. 


2001) 


(2002 


) (red dotted line) 



lifetime for three dif f erent weak lensing estimates of the Coma cluster mass ( IKubo et al.ll2007 



Gavazzi et al.ll2009t lOkabe et al.l 120101 ) while the green dashed line shows constraints if the 



mass is estimated based on the velocity dispersion of the cluster galaxies (ILokas fc Mamon 

iooi ). 
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Fig. 3. — Constraints on the decay lifetime as a function of mass from dwarf spheroidal 
galaxies for a bb final state (left panel) and a fi + n~ final state (right panel) based on the 
95% confidence level upper limits on the gamma-ray flux from 11 months of Fermi-L AT 
observations. Here we conservatively assume no IC gamma-ray emission. 
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Fig. 4. — Constraints on the dark matter decay lifetime as a function of particle mass for 
the Ursa Minor dwarf spheroidal galaxy for an assumed final state. Left: Constraints 

including IC emission for different choices of the diffusion coefficient Dq. The black line 
conservatively shows the constraints for no IC emission, the blue lines include IC emission 
for two diffusion coefficients (Do = 10 28 cm 2 s : dot-dashed line and Do = 10 29 cm 2 s : 
dashed line) bounding those d erived for th e Milky Way and similar to what was assumed 
for dark matter annihilation in lAbdo et al.l (j2010bl ). and the red line shows the constraints 
in the limit of no diffusion. Right: Constraints for a diffusion coefficient of D = 10 28 cm 2 
s _1 for different assumed masses for Ursa Minor. The solid line reproduces the constraints 
using the conservative assumption of a mass corresponding to the half-light radiu s, M\ii. For 



compar ison, we also show the constraints using the NFW profile parameters from I Ab do et al. 



( I2010bl ) based on fits to the Ursa Minor stellar velocity data and calculating Jd within a radius 
of 1 degree (dot-dashed line). The dashed lines show constraints on the decay lifetime for an 
order of magnitude variation in the assumed mass; these masses bracket the expected mass 
within a radius of 1 degree for dark mater density profiles fitting the stellar data for Milky 
Way dwarf sphero idal galaxies which are all consistent with having the same total halo mass 
( jWolf et al.ll2010l . see for example their Figure 3). 
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Fig. 5. — Constraints on the dark matter decay lifetime as a function of particle mass for 
the M31 galaxy, for a bb final state (left) and for a fi + n~ final state (right). In the left 
panel, we show with a black line the constraints obtained employing a NFW profile with 
parameters from iGeehan et al.l (120061 ) (assuming a 1 degree projected radius) and the 2-er 
upper limit to the gamma-ray emission with an assumed point-like source spatial distribution; 
th e green dashed line uses the same gamma-ray limit, but assumes the iso-thermal profile 
of iKerins et al.l ( 120011 ). Finally, for the blue dotted line we employ the M31 virial mass and 
the Fermi 2-cr upper limit to the gamma-ray emission with an assumed extended source 
distribution. In the right panel, we show the constraints on a final state, for diffusion 

coefficients set to D = 10 28 cm 2 /s (solid lines) and D = 10 29 cm 2 /s (dashed lines); the 
black lines refer to the NFW plus point-like source assumptions, while the red lines to the 
case where we employ the virial mass and the extended source gamma-ray flux. 



